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,H ' Abstract 

& . We have used the Japanese VLBI array VERA to perform high-precision astrometry of an H2O maser 

source in the Galactic star-forming region NGC 281 West, which has been considered to be part of a 300-pc 
I I superbubble. We successfully detected a trigonometric parallax of 0.355±0.030 mas, corresponding to a 

j> ■ source distance of 2.82±0.24 kpc. Our direct distance determination of NGC 281 has resolved the large 

' distance discrepancy between previous photometric and kinematic studies; likely NGC 281 is in the far 

C*~) . side of the Perseus spiral arm. The source distance as well as the absolute proper motions were used to 

demonstrate the 3D structure and expansion of the NGC 281 superbubble, ~650 pc in size parallel to the 
Galactic disk and with a shape slightly elongated along the disk or spherical, but not vertically elongated, 
\& • indicating the superbubble expansion may be confined to the disk. We estimate the expansion velocity of 

the superbubble as ~ 20 km s _1 both perpendicular to and parallel to the Galactic disk with a consistent 
25 , timescale of ~ 20 Myr. 

Key words: Galaxy: kinematics and dynamics — ISM: bubbles — ISM: Hll regions — ISM: individual 
(NGC 281) — masers (H 2 0) 
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JH I 1. Introduction pershclls, many of which were identified earlier by Hcilcs 

(1979, 1984) based on Hi 21-cm line surveys by Weaver 

A paradigm shift has been occurring in recent years in and Williams (1973) and Heiles and Habing (1974). Kiss, 

the observational study of the interstellar medium (ISM) Moor, and Toth (2004) later identified 145 large-scale 

in the Galaxy and in external galaxies, where the ISM is loop- or arc-like structures in the 2nd Galactic quadrant 

now recognized as being of fundamental importance not (i.e. for galactic longitudes 90° < I < 180°) on IRAS-b&scd 

only in the life cycle of stars but also in galactic evolution, far- infrared maps. Recent results for the 1st, 3rd, and 4th 

The present understanding of the ISM draws a dynamic Galactic quadrants by Konyves et al. (2007) increased the 

and violent picture of the gas in the vast interstellar space, number of identified loops to 462. Similar objects are also 

rather than a static one envisioned previously. Recent ob- found in nearby galaxies, including the Local Group spi- 

servations with higher angular resolution have revealed rals Andromeda (M31; Brinks & Bajaja 1986) and M33 

the existence of dynamic and diverse gas structures in the (Deul & den Hartog 1990), Local Group irregulars includ- 

Galaxy that cannot be explained by the traditional evolu- ing the Large Magellanic Cloud (LMC; Chu et al. 1993, 

tionary cycle of the ISM, requiring new insights toward a 1995; Kim et al. 1998, 1999), the Small Magellanic Cloud 

comprehensive understanding of the ISM (e.g., Reynolds (SMC; Staveley-Smith et al. 1997; Stanimirovic et al. 

1997, 2002; Dickey 2001). Extreme activity of the ISM is 1999) and IC 10 (Shostak & Skillman 1989), large spi- 

indicated by large-scale circular or arc-like structures (up ral galaxies such as M101 (Kamphuis et al. 1991) and 

to a few kiloparsecs in scale) called superbubbles or su- NGC 891 (Howk & Savage 1997, 2000), and dwarf galax- 
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Fig. 1. A wide-field optical image of the star-forming region NGC 281 taken with the WIYN 0.9-meter at Kitt Peak National 
Observatory (courtesy of T.A. Rector/University of Alaska Anchorage and WIYN/AURA/NSF). The orientation of the image is 
north up and east left. The red nebulosity of the Hn region NGC 281 is obscured in its southwestern quadrant by the adjoining 
molecular cloud NGC 281 West, where the H2O maser source is located. 



ies such as M82 (Matsushita et al. 2000; Matsumoto et al. 
2001), IC 2574, DDO 47 (Walter & Brinks 1999, 2001) 
and Holmberg II (HoII; Puchc et al. 1992). 

The most widely accepted theory for the origin of su- 
perbubbles suggests that supernova explosions and strong 
stellar winds from OB associations in a galactic disk blow 
the surrounding gas up into the halo, carving out exten- 
sive cavities in the ISM (see Tenorio-Tagle & Bodenheimer 
1988). Superbubbles are therefore believed to play a vital 
role in disk-halo interactions, serving as a conduit for mat- 
ter and energy from the disk into the halo in the form of 
galactic 'fountains' (Shapiro & Field 1976) or 'chimneys' 
(Norman & Ikeuchi 1989). 

In the Galaxy, only a small number of chimneys (or also 
called 'worms') or fragmenting supcrbubbles/supcrshells 
(on several hundred to one thousand parsec scales) are 
known and well-studied: e.g., the Orion-Eridanus su- 
pcrbubble (Cowie et al. 1979; Reynolds & Ogden 1979; 
Boumis et al. 2001; Welsh et al. 2005), the Cygnus 
superbubble (Cash et al. 1980), the Stockert chimney 
(Muller et al. 1987; Kundt & Muller 1987; Forbes 2000), 
the Aquila supershell (Maciejewski et al. 1996), the 
Scutum supershell (Callaway et al. 2000; Savage et al. 
2001), GSH 277+00+36 (McClure-Griffiths et al. 2000), 
GSH 242-03+37 (McClure-Griffiths et al. 2006), the W4 



chimney /superbubble (Normandeau et al. 1996; Reynolds 
et al. 2001; Terebey et al. 2003; Madsen et al. 2006; West 
et al. 2007), the Ophiuchus superbubble (Pidopryhora 
et al. 2007), and the NGC 281 superbubble (Megeath et al. 
2002, 2003; Sato et al. 2007). Superbubbles are also a key 
target of research for understanding their impact on the 
history of star formation in the galactic disk, including 
but not limited to supcrbubble-triggered star formation 
(e.g., NGC 281 as explained later in the section). 

In either the Galaxy or external galaxies, very little is 
known about the dynamics and energetics of superbub- 
bles/chimneys or about their influence on the formation 
of the halo or on the star formation history in the galactic 
disk. A difficulty in revealing the structure and kinemat- 
ics of superbubbles stems from the limited information ob- 
tainable from two dimensions of sky positions and velocity 
information only along the line of sight. A direct observa- 
tion of the three-dimensional (3D) motion and structure 
of a superbubble has not been easy due to the high astro- 
metric accuracy required. However, Sato et al. (2007) re- 
ported the first VLBI (Very Long Baseline Intcrferometry) 
observations that obtained direct evidence for expanding 
motions of a superbubble by measuring the proper mo- 
tions. 

In Sato et al. (2007), we observed the H2O maser emis- 



Distance to NGC 281 



3 



sion in NGC 281 West over 6 months with VERA (VLBI 
Exploration of Radio Astrometry), a new Japanese radio 
telescope array dedicated to phase-referencing VLBI as- 
trometry (e.g., Honma et al. 2000, 2005, 2007; Kobayashi 
et al. 2003; Imai et al. 2006, 2007; Hirota et al. 2007, 
2008a, 2008b; Nakagawaet al. 2008; Kim et al. 2008; Choi 
et al. 2008). These observations revealed the systemic mo- 
tion of the NGC 281 superbubble away from the Galactic 
plane at a velocity of 20—30 km s^ 1 . At its estimated 
heliocentric distance of 2.2—3.5 kpc, NGC 281 most likely 
originated in the Galactic plane and has been blown out 
by a superbubble expansion (see §1.1 and Sato et al. 2007). 

A precise and direct determination of the distance to 
NGC 281 is desirable and essential for a better under- 
standing of the region, including further description of 
the 3D structure and motion of the superbubble to reveal 
the origin, the energetics and the timescale of the super- 
bubble. These are the primary reasons why we decided 
to complete our VERA observations to directly measure 
the parallactic distance to NGC 281 via an H2O maser 
source. The astrometric determination of the distance 
to NGC 281 is also important for studying early- type 
(O— B) high-mass (M> IOMq) stars and stellar evolution 
in combination and in comparison with photometric stud- 
ies of the OB stars in NGC 281 (e.g., Henning et al. 1994; 
Guetter & Turner 1997). In this paper, we present the 
results of our parallax measurements with VERA over 18 
months, which enables us to further discuss the 3D struc- 
ture, expanding motion, the origin and the timescale of 
the NGC 281 superbubble. 

1.1. NGC 281 Superbubble 

The NGC 281 region (a 200 o =00 h 52 m , c5 200 o = +56°34' 
or / =123°. 07, b = — 6.°31) provides an excellent laboratory 
for studying in detail the cycle of the ISM and its impact 
on both star formation and Galactic disk-halo interaction 
through superbubbles. At its estimated heliocentric dis- 
tance of 3 kpc, this region is remarkably located far above 
(~ 300 pc) the midplane of the Perseus arm of the Galaxy. 
Of special interest in the NGC 281 region is a possibility of 
triggered star formation occurring on two different scales: 
large-scale (~300 pc) supperbubblc-triggercd formation of 
the first OB stars (Megeath et al. 2002, 2003), and the se- 
quential and ongoing triggered star formation (~1 — 10 pc) 
in an adjoining giant molecular cloud (NGC 281 West) 
through interaction with an Hll region (the NGC 281 neb- 
ula) excited by the first OB stars (Elmegreen & Lada 1978; 
Megeath & Wilson 1997). 

Figure 1 is an optical image of the region, which clearly 
shows the red nebulosity of the Hll region NGC 281 (also 
known as Sharplcss 184) of diameter 20 arcmin (~ 17 pc), 
in which an early- type cluster called IC 1590 is embed- 
ded (Guetter & Turner 1997). The brightest member of 
the cluster IC 1590 is an O-type Trapezium-like system 
HD 5005 (or also called ADS 719), whose component stars 
HD 5005ab (unresolved), HD 5005c, and HD 5005d have 
spectral types of 06.5 V, 08 V, and 09 V, respectively 
(Walborn 1973; Abt 1986; Guetter & Turner 1997). 

The southwestern quadrant of the NGC 281 nebula is 



obscured (as seen in figure 1) by the adjoining molecular 
cloud NGC 281 West in front of the nebula. Ongoing star 
formation in the NGC 281 West cloud is indicated by the 
H2O maser emission and IRAS sources within this cloud 
near the clumpy interface between the H II region and the 
cloud. This generation of stars may have been triggered 
by interactions with the Hll region (Elmegreen & Lada 
1978; Megeath & Wilson 1997). The NGC 281 molec- 
ular cloud complex that surrounds the Hll region (in- 
cluding the NGC 281 West and East clouds) was mapped 
in CO emission lines by Lee and Jung (2003). The cen- 
tral radial velocity of the NGC 281 West molecular cloud, 
Vlsr = 31 km s _1 (Lee and Jung 2003), agrees well with 
that of the H2O maser emission in the cloud (Sato et al. 
2007). On a larger scale, Megeath et al. (2002, 2003) iden- 
tified the CO molecular cloud complex as being on an Hi 
loop extending over 300 pc away from the Galactic plane 
toward decreasing galactic latitude. In an I vs. Vlsr dia- 
gram (i.e. a plot of galactic longitude vs. observed radial 
velocity), the CO molecular clouds appear to be part of 
a broken ring of diameter 270 pc, expanding at a veloc- 
ity of 22 km s~ J parallel to the Galactic plane (Megeath 
et al. 2002, 2003). Megeath et al. (2002, 2003) suggested 
that these clouds were formed in a fragmenting superbub- 
ble, which triggered the formation of the first OB stars, 
and that these OB stars then ionized the surrounding gas 
which subsequently triggered ongoing star formation in 
the neighboring cloud. 

In Sato et al. (2007), we derived the absolute proper 
motions of the H2O maser features in NGC 281 West by 
phase-referencing VLBI observations with VERA over 6 
months. Wc found 10 maser features, highly likely excited 
by two spatially distinct young stellar objects (YSOs), 
are systematically moving toward southwest, in agree- 
ment with the direction of the proper motion of HD 5005 
in the Hll region measured with larger uncertainties by 
Hipparcos (Perryman et al. 1995, 1997). Our results with 
VERA revealed the systemic motion of the NGC 281 re- 
gion away from the Galactic plane, most likely having 
originated in the Galactic plane and been blown out by 
the expansion of the NGC 281 fragmenting superbubble. 

The uncertainties in our velocity estimates are largely 
due to the large distance uncertainty of NGC 281, despite 
many measurements, e.g., photometric distances of 2 kpc 
(Sharpless 1952, 1954), 3.68 kpc (Cruz-Gonzalez et al. 
1974), 3.5 kpc (Henning et al. 1994) and 2.94±0.15 kpc 
(Guetter & Turner 1997) and kinematic distances of 
2.2 kpc (Gcorgelin & Georgelin 1976; Roger & Pedlar 
1981) and 3 kpc (Lee & Jung 2003). In Sato et al. (2007), 
we adopted a value of 2.9 kpc as the distance to NGC 281, 
based on a recent estimate of 2.94±0.15 kpc derived by 
Guetter and Turner (1997) through photometry of 279 
individual stars in and about the young cluster IC 1590 
embedded in the NGC 281 nebula. 

The most reliable and direct way to determine the dis- 
tance (d) of an astronomical object is to measure the 
trigonometric parallax (71") of the object due to the Earth's 
motion around the solar system barycenter, where the dis- 
tance is given by rf(pc) = 1/tt("). The astrometric deter- 
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mination of the distance to NGC 281 provides valuable 
information for studying the formation and evolution of 
early-type (O— B) high-mass (M > IOMq) stars in com- 
parison with photometric studies and distances, for exam- 
ple, by Guetter & Turner (1997). 

In the present study, we report on our successful deter- 
mination of the parallactic distance to NGC 281, which 
reveals the 3D structure of the NGC 281 fragmenting su- 
perbubble, and investigate the origin and timescale of the 
formation of such a large-scale (~650 pc) structure. 

2. VERA Observations and Data Reduction 

The VERA observations were carried out over 13 epochs 
spaced approximately at monthly to bimonthly intervals. 
The first 6 epochs reported in Sato et al. (2007) were on 
2006 May 14, July 21, August 3, September 5, October 
25, and November 18 (days of year 134, 202, 215, 248, 
298, and 322, respectively). Seven new observations were 
conducted on 2006 December 11, 2007 January 22, March 
19, May 1, July 27, September 30, and October 28 (days of 
year 345, 387, 443, 486, 573, 638, 666 as counted from 2006 
January 0). Note that the dates stated here are universal 
time (UT) start times, while for the data analyses in the 
following section we use more precisely the median hour 
of each observation, which lasted 7—9 hours. 

The observational procedures are as described in Sato 
et al. (2007) and as follows. We observed two sources si- 
multaneously in the dual-beam mode of VERA for phase 
referencing (e.g., Honma et al. 2003), and the real-time in- 
strumental phase difference data between the two beams 
were taken with an artificial noise source in each beam and 
recorded for calibration (Kawaguchi ct al. 2000; Honma 
ct al. 2008a). We simultaneously observed the H2O 
maser source in the NGC 281 West cloud with one of the 
two extragalactic position-reference quasars J0047+5657 
(0°.84 separation; at a 20 oo =00 h 47 m 00. s 428805, <5 200 o = 
+56°57'42."39479) and J0042+5708 (1°.50 separation; 
at a 20 oo =00 h 42 m 19. s 451727, d 20 oo = +57°08'36."58602). 
Positions are from Bcasley et al. (2002). We alternated be- 
tween the two quasars typically every 10 minutes. We ob- 
served a bright calibrator source hourly: J0319+4130 (= 
3C 84) at the first epoch (2006 May 14) and J2232+1143 
(= CTA 102) for the other 12 epochs. The data correlation 
was performed with the Mitaka FX correlator (Chikada 
et al. 1991) with frequency and velocity resolutions for 
the maser lines of 15.625 kHz and 0.21 km s _1 , respec- 
tively. The observed frequencies of the maser lines were 
converted to radial (line-of-sight) velocities with respect 
to the local standard of rest (LSR), Vlsr, using a rest fre- 
quency of 22.235080 GHz for the H 2 6i6-5 2 3 transition. 

Visibility calibration and imaging were performed in 
a standard manner with the NRAO Astronomical Image 
Processing System (AIPS) package. The visibilities of all 
H 2 maser channels at each epoch were phase-referenced 
to each of the position-reference sources J0047+5657 and 
J0042+5708 independently by applying the phase solu- 
tions from fringe fitting (AIPS taks FRING) for each 
position- reference source to the H 2 maser channels for 



the corresponding observation time and frequencies. The 
instrumental phase difference between the two beams was 
also corrected using the recorded phase-calibration data 
described above. We also calibrated the error in the vis- 
ibility phase caused by the Earth's atmosphere based on 
GPS measurements of the atmospheric zenith delay due 
to tropospheric water vapor (Honma et al. 2008b). After 
calibration, we made spectral-line image cubes using the 
AIPS task IMAGR around the maser spots with 512 pix- 
els x 512 pixels of size 0.05 mas. The synthesized beam 
had an FWHM beam size of 1.3 mas x 0.8 mas with a po- 
sition angle of —43° East of North. Least-squares fittings 
were performed to simultaneously solve for the sinusoidal 
parallax curve and linear proper motion in the directions 
of right ascension (RA) and declination (Dec). 

3. Results 

3.1. Parallax Measurements 

Figure 2 shows the position measurements of the 
brightest maser spot at Vlsr = —32.1 km s _1 (feature 
#4 of 10 features identified in Sato et al. 2007) using 
J0047+5657 (blue) and J0042+5708 (red) as position- 
reference sources independently. Figures 2a and 2b show 
the positional offsets of the maser spot in the eastward 
(X, i.e. (RA) cos(Dcc)) and northward (Y, i.e. Dec) di- 
rections, respectively, relative to the reference position at 
the origin at a 200 o =00 h 52 m 24. s 70081 for X and 5 20 oo = 
+56°33'50."5274 for Y, as a function of time (day of year, 
as counted from 2006 January 0). The constant offsets of 
~ 0.80 mas seen in each of X and Y between the two 
measurements using different position-reference sources 
J0047+5657 and J0042+5708 arc due to the uncertain- 
ties in the absolute positions of these reference sources 
(0.64 mas for J0047+5657 and 0.82 mas for J0042+5708; 
Beasley et al. 2002) and are solved for in the multi-epoch 
parallax and proper- motion fittings. The best-fit models 
are also plotted with gray lines for proper motions and 
black curves for the sum of sinusoidal parallax and the 
linear proper motions. Table 1 summarizes the fitted dis- 
tances and proper motions. Figures 2c is a plot of the 
spot position in Y vs. X and its movement during the 
observation epochs. 

The obtained values of the parallax in X are: tt = 
0.346 ± 0.052 mas (corresponding to a distance of d = 
2.89 ± 0.43 kpc) using J0047+5657 and tt = 0.334 ± 
0.047 mas (corresponding to a distance of d = 2.99 ± 
0.42 kpc) using J0042+5708. Combining these indepen- 
dent measurements using two different position-reference 
sources, we have obtained as a final value in X (RA) di- 
rection: tt = 0.340 ± 0.034 mas (d = 2.94 ± 0.30 kpc). 

Similarly, we have obtained the parallax in Y to be: tt = 
0.400 ± 0.095 mas (d = 2.50 ± 0.59 kpc) using J0047+5657 
and tt = 0.412 ± 0.092 mas (d = 2.43 ± 0.54 kpc) using 
J0042+5708. Combining these two measurements yields 
the final value in Y (Dec) direction to be: tt = 0.406 ± 
0.065 mas (d = 2.46 ± 0.39 kpc). Here we have estimated 
the associated uncertainties from the uniformly-weighted 
standard deviation a from the least-squares fitting in X 
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Fig. 2. Position measurements of the peak maser spot of feature 4; using J0047+5657 (blue) and J0042+5708 (red) as position-ref- 
erence sources. The coordinate origin is at a 2 ooo =00 h 52 m 24. s 70081, <5 2 ooo = +56°33'50f 5274. (a) and (b) show the eastward 
(X =RA cos(Dec)) and northward (Y =Dec) positional offsets in mas, versus time (day of year, as counted from 2006 January 0), 
respectively. Best-fit models are also plotted with gray lines for proper motions and black curves for the sum of sinusoidal parallax 
and the linear proper motions, (c) The trajectory of the positional offsets X versus Y (in mas) on the sky. Best-fit models are 
also plotted, with cyan crosses showing the expected positions on each day of observation. The associated error bars in each panel 
indicate the standard deviation in X or Y from the least-squares fits (see text). 
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Table 1. Fitted distances and absolute proper motions for the peak maser spot of feature 4. 



X , Y Reference ir [mas] d [kpc] a [mas] fix [mas yr r ] /iy [mas yr x ] 

(1) (2) (3) (4) (5) (6) (7) 

X J0047+5657 0.346 (0.052) 2.89 (0.43) 0.113 -2.60 (0.07) 

(RA) J0042+5708 0.334 (0.047) 2.99 (0.42) 0.104 -2.65 (0.07) 

Combined 0.340 (0.034) 2.94 (0.30) 0.106 -2.63 (0.05) 

Y J0047+5657 0.400 (0.095) 2.50 (0.59) 0.188 -1.97 (0.12) 

(Dec) J0042+5708 0.412 (0.092) 2.43 (0.54) 0.183 -1.75 (0.11) 

Combined 0.406 (0.065) 2.46 (0.39) 0.181 -1.86 (0.08) 

X fc Y Combined 0.355 (0.030) 2 82 (0.24) 0.148 

(1) Directions of position measurement: eastward (X) and northward (Y). Data were uniformly weighted for each 
epoch. For the X & Y direction, data in the two directions X and Y were error-weighted with the standard deviation 
in the fitting in each direction. (2) Position-reference sources used for the measurement. 'Combined' indicates the 
simultaneous fitting of two independent position-reference results using J0047+5657 and J0042+5708. (3) Measured 
parallaxes (in mas) with associated uncertainties in parentheses. (4) Distances corresponding to the parallax (in 
kpc) with associated uncertainties in parentheses. (5) The standard deviation of the post-fit residuals (in mas) from 
the least-squares fit. (6)(7) Absolute proper motions in X and Y (in mas yr - 1 ) with associated uncertainties in 
parentheses. 



and Y, which were ax = 0.104 — 0.113 mas and ay = 
0.181 — 0.188 mas, respectively (see table 1 for the individ- 
ual value of the standard deviation in each measurement). 
Random errors in individual position measurements esti- 
mated with the AIPS Gaussian-fit task JMFIT caused by 
the noise (which are approximated by the beam HWHM 
divided by the signal-to-noise ratio) were ~ 0.010 mas, 
which is thus not the main cause of the deviations ax 
and ay from the fits. We will discuss the possible error 
sources in §4.1 in more detail. However, those errors are 
difficult to measure quantitatively, and we therefore esti- 
mated the astrometric errors by the standard deviations 
of the residuals from the fits and plotted the standard 
deviations as the error bars of each point in figure 2. 

Weighting by ax and ay , we simultaneously fitted the 
position measurements of the maser spot in both X and 
Y altogether, including all data using J0047+5657 and 
J0042+5708, to obtain the final result for the parallax 
measurements to be: ir = 0.355 ±0.030 mas, which corre- 
sponds to a distance to NGC 281 of d = 2.82 ±0.24 kpc. 
We have therefore successfully determined the parallac- 
tic distance to NGC 281 as far as ~ 3 kpc with a preci- 
sion better than 10%. Our distance estimate of NGC 281, 
d = 2.82 ± 0.24 kpc, agrees well with the recent photo- 
metric distance of 2.94±0.15 kpc by Guetter and Turner 
(1997), but excludes the kinematic distance of 2.2 kpc by 
Georgelin and Georgelin (1976) and by Roger and Pedlar 
(1981) and most distance estimates in earlier studies. 

3.2. Absolute Proper Motions 

As listed in table 1, we fitted the position measure- 
ments of 13 epochs simultaneously for both the parallax 
7T and proper motions (fix, My) (in RA and Dec, respec- 
tively) of the peak maser spot of feature 4 (the bright- 
est feature in the maser component C3, which is one of 
the two spatial components CI and C3 with maser emis- 
sion detected in Sato et al. 2007). For the proper mo- 
tions, we obtained (-2.60±0.07 -1.97±0.11) mas yr" 1 us- 
ing J0047±5657, and (-2.65±0.07, -1.74±0.11) mas yr" 1 



using J0042+5708. The error-weighted mean proper 
motions of this spot (which are the same as un- 
weighted mean in this case of equal errors) are: 
{fix, My)=(-2.63±0.05, -1.86±0.08) mas yr" 1 . We con- 
sider this mean as revised values (called 'modified' here- 
after) of the error-weighted mean proper motions (denoted 
by 4w) of feature 4 reported in Sato et al. (2007). 

Table 2 shows a comparison of the absolute proper mo- 
tions derived for feature 4 between the revised values in 
the present study and the previous values in Sato et al. 
(2007). In Sato et al. (2007), we derived absolute proper 
motions of feature 4 (in C3) and features 9 and 10 (in 
CI, the other one of the two spatial maser components; 
see Sato et al. 2007) assuming the distance to NGC 281 
to be d = 2.9 kpc (adopting the photometric distance by 
Guetter and Turner 1997) and subtracting the expected 
parallax in the position measurements during the first 6 
observation epochs. 

The revised values of fix for feature 4 in the present 
study agree with the previous values in Sato et al. (2007) 
within the margin of errors {la), however, for fiy we find 
deviations of the revised values from the previous ones 
larger than the errors (> 2a) estimated in Sato et al. 
(2007). This means in Sato et al. (2007) we overestimated 
/iy of feature 4 and underestimated the errors in the lin- 
ear least-squares fittings (after subtracting the assumed 
parallax) . As seen in figure 2 based on the present longer- 
time observations, the overestimate of the previous fiy 
values for feature 4 is due to a large position deviation 
from the best-fit at the first (±2<r) epochs, which yielded 
a larger slope. The underestimate of the errors occurcd 
mostly due to non-gaussian errors larger than 0.1 mas 
(larger than thermal errors by an order of magnitude) 
in the position measurements, owing to astrometric error 
sources that are discussed in more detail in §4.1. However, 
as will be described in §4.2, the scientific conclusions re- 
main essentially unchanged after the modifications to the 
proper motions. 

By adding the mean relative motion of all 8 fea- 
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Table 2. Comparison of the absolute proper motions of feature 4 between the present study and Sato et al. (2007). 



Reference 
(1) 


fix [mas yr 1 ] 
Present study 
(2) 


fix [mas yr L ] 
Sato et al. (2007) 
(3) 


fly [mas yr l ] 
Present study 
(4) 


fly [mas yr L ] 
Sato et al. (2007) 
(5) 


J0047+5657 
J0042+5708 


-2.60 (0.07) 
-2.65 (0.07) 


-2.87 (0.26) 
-2.92 (0.23) 


-1.97 (0.12) 
-1.75 (0.11) 


-2.78 (0.37) 
-2.54 (0.35) 


Combined 


-2.63 (0.05) 


-2.89 (0.18) 


-1.86 (0.08) 


-2.65 (0.26) 



(1) Position-reference sources used for the measurement. 'Combined' indicates the simultaneous fitting of two inde- 
pendent position-reference results using J0047+5657 and J0042+5708. (2) (3) Absolute proper motions of feature 4 
derived in the present study and in Sato et al. (2007), respectively, in X direction (eastward) (in mas yr - 1 ) with 
associated uncertainties in parentheses. (4)(5) Same as (2)(3) but in Y direction (northward). 



tures (including feature 4) in C3 with respect to the 
reference feature 4, (fl x , ^)=(0.46±0.33, 0.42±0.31) 
mas yr -1 as given in Sato et al. (2007), to this mod- 
ified absolute motion of feature 4, the mean absolute 
proper motion of these features in C3 (denoted by 
C3m in Sato et al. 2007) is then also revised to be: 
(fix, /iy) C 3 =(-2.17±0.33, -1.44±0.32) mas yr" 1 . 

Our new distance estimate of NGC 281 by parallax 
measurements, d = 2.82 kpc, also introduces a slight 
modification to the measured absolute proper motions 
of features 9 and 10 on the same data of 6 epochs 
(as these features were not persistent through the to- 
tal 13-cpoch observations with high enough signal-to- 
noise ratios for the parallax and proper motion fit- 
tings). The modified absolute proper motions (fix, fly) 
of feature 9 are (-4.85±0.11, -1.73±0.46) mas yr" 1 
using J0047+5657, and (-4.74±0.28, -1.44±0.41) 
mas yr -1 using J0042+5708. Similarly we obtained 
the absolute proper motion (fix, fiy) of feature 10 
to be: (-2.00±0.23, -4.06±0.79) mas yr" 1 using 
J0047+5657, and (-1.89±0.38, -4.14±0.69) masyr" 1 us- 
ing J0042+5708. Here the associated uncertainties in fix 
and fly were estimated from the standard deviations from 
the linear least-squares fits in X and Y positional offsets 
(after the parallax for d = 2.82 kpc is subtracted), respec- 
tively. The error- weighted mean proper motions of each of 
these features in CI (denoted by 9w and lOw in Sato et al. 
2007) arc [fix, £ty)=(-4.84±0.10, -1.57±0.31) mas yr" 1 
for feature 9 and (fi x , fiy)=(- 1-97 ±0.20, -4.11±0.52) 
mas yr -1 for feature 10. Then the unweighted mean of 
the absolute motions of these two features in CI (denoted 
by Clm in Sato et al. 2007) is also slightly modified to be 
(fix, ffvhi =(-3.41, -2.84) mas yr" 1 . 

4. Discussion 

4-.1. Astrometric Error Sources 

In the present work, we measured the annual parallax of 
NGC 281 West to be 0.355±0.030 mas, which corresponds 
to a distance of 2.82 ± 0.24 kpc. In this section, we first 
discuss possible error sources in our parallax and proper- 
motion measurements. 

As already mentioned in §3.1, the true uncertainties in 
the measurements were estimated from standard devia- 
tions of the post-fit residual from the least-squares fits 
(as listed in table 1), because thermal errors due to noise 



in individual position measurements (~ 0.010 mas) are 
far smaller than the deviations from the fits. The likely 
sources for such deviations are difficult to measure quan- 
titatively, which led to our error underestimate in the de- 
rived absolute proper motions of feature 4 in Sato et al. 
(2007) as described in §3.2. The reference sources did 
not show any resolved structure such as additional com- 
ponents or structural change, and we obtained similar re- 
sults. Therefore we do not consider the contribution of 
the reference source as a dominant error source. 

One of the likely causes for position uncertainties of 
ax = 0.104-0.113 mas and a Y =0.181-0.188 mas (see 
table 1) is mis-modeling of troposphcric zenith delay, as 
also reported in previous astrometric measurements with 
VERA (e.g., Hirota et al. 2007; Nakagawa et al. 2008; 
Honma et al. 2008b). The water vapor in the tropo- 
sphere introduces optical path differences through the at- 
mosphere between the target maser source and reference 
sources because of the elevation angle difference between 
the sources. Therefore, even after the phase referencing 
calibration, there remain residuals of troposphcric zenith 
delay due to the water vapor, which are difficult to mea- 
sure precisely and are one of the most serious error sources 
in the VLBI astrometry in the 22-GHz band. Nakagawa 
et al. (2008) and Honma, Tamura & Reid (2008b) give 
detailed discussion on astrometric errors caused by such 
residuals of troposphcric zenith delay calibration with 
VERA and report probable errors, ax and ay, of 0.05- 
mas and 0.10-mas levels, respectively, for source declina- 
tion of 60°. Our parallax determination with VERA is 
no exception in having astrometric errors larger in decli- 
nation (Y) than in right ascension (X) (i.e. ay > ax) as 
the same tendencies reported in other observations with 
VERA (e.g., Hirota et al. 2007). Larger errors in Y (Dec) 
than in X (RA) arc explained by the tendency of tro- 
posphcric zenith delay errors to have more severe effects 
in the Y direction (although depending of the separation 
angle between the target and reference sources and the 
declination of the sources) and are as expected by the 
simulation in Honma, Tamura and Reid (2008b). 

However, we do not consider the residuals of tropo- 
sphcric zenith delay as the dominant cause of errors in our 
parallax measurements. As the simulations by Honma, 
Tamura and Reid (2008b) clearly show, the astrometric 
errors due to the zenith delay residuals for the NGC 281 
case at its high declination (^ooo — +56°. 5) and with 
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separation angle (SA) smaller than 1° for J0047+5206 
(SA=0°.84) are expected to be as good as 20 — 30 fias 
(i.e. 0.020 — 0.030 mas) in each of X and Y directions, 
which do not account for all the standard deviations 
cr x = 0.104-0.113 mas and o Y = 0.181 - 0.188 mas in 
our measurements. Also, as figure 2 clearly demonstrates, 
using two different reference quasars with different sepa- 
ration angles (SA=0°.84 for J0047+5206 and SA=1°.50 
for J0042+5708) did not result in remarkable differences 
in the astrometric errors for the two cases. This certainly 
excludes the simple zenith delay error effect as the main 
source of our astrometric errors. 

The most likely major source of errors in the present po- 
sition measurements is a variation of the intrinsic maser 
structure. Most maser spots are unresolved in VLBI ob- 
servations, however evidence has been shown for intrinsic 
maser structure smaller than 1 AU on submilliarcsecond 
scales in previous maser studies (e.g., Fish & Sjouwerman 
2007). A maser spot likely consists of both unresolved 
substructure and observed larger structure which arc in- 
distinguishable (e.g., Fish & Sjouwerman 2007), and thus 
a variation in maser substructure can lead to fluctuation of 
the peak position, limiting the accuracy in position mea- 
surements of a maser spot. This effect of the maser spatial 
structure is also seen for other observations in parallax 
measurements with up to f0% errors, independently of 
the source distance, for example, in Orion KL by Hirota 
et al. (2007). Errors due to the maser structure can be 
reduced by using multiple maser features for derivation of 
a parallax. In our case of NGC 281 West, however, only 
one bright maser feature (above 10 Jy beam -1 ) that was 
used in the present study was found with a good signal- 
to-noise ratio for the precise parallax measurements, even 
persistent for the all observation epochs. 

According to discussions by Honma et al. (2007), other 
possible astrometric errors in VERA observations that 
arise from the uncertainties in the antenna-station posi- 
tions, delay models, and path length differences due to 
ionosphere are smaller by at least an order of magnitude 
than the uncertainties due to troposphcric zenith delay 
residuals. Therefore, we conclude that the main error 
sources of our astrometric measurements is the variabil- 
ity of maser feature structure, and the contribution by 
troposphcric zenith delay residuals is also nonnegligiblc. 

4-2. NGC 281 Superbubble vs. Galactic Rotation 

Using our new results from parallax and proper-motion 
measurements of the H2O maser source in NGC 281 West, 
we now revise the motion of the NGC 281 region with 
respect to the Galactic rotation as reported in Sato et al. 
(2007). 

In an analogous manner to Sato et al. (2007), we show 
the modified systemic motion of the NGC 281 region away 
from the Galactic plane, traced by our H2O maser obser- 
vations, as a red arrow in figure 3a, relative to a frame 
rotating with the Galactic rotation. Figures 3a, b and c 
are an image of the NGC 281 superbubble plotted as b 
vs. I (velocity-integrated for the Perseus-arm line-of-sight 
velocity range, Vlsr = —60 to —25 km s _1 ) and Vlsr vs. 



I (latitude-integrated for the galactic latitude range in- 
dicated by a cyan rectangle in figure 3a) diagrams after 
Megeath et al. (2003) and based on Hi 21-cm line data 
by Hartmann and Burton (1973) and 12 CO(J= 1 — 0) line 
data by Dame, Hartmann, and Thaddeus (2001). At the 
position of NGC 281 West, the Galactic plane lies almost 
parallel to the east- west (RA) direction with a position an- 
gle of 90. °2, so that we can regard the motions in X (RA) 
and Y (Dec) directions as those parallel (I) and perpen- 
dicular (6; with only small projection effect at 6 = — 6.°31) 
to the Galactic plane, respectively. 

Subtracting the expected apparent motions of NGC 281 
due to the (simulated) Galactic rotation (with re- 
spect to the LSR) and the solar motion effect, 
(/XGX; Mgf)=(— 3.27, —0.29) mas yr _1 (including the ap- 
parent motion arising from the nonzero galactic latitude 
of NGC 281; 6 = -6.°31) and (fi sx , msy) = (0.84, -0.53) 
mas yr _1 , respectively, from the measured mean proper 
motion of the H2O maser source in NGC 281 West yields 
the systemic motion of the region with respect to the ro- 
tation of the Galaxy. Here we have adopted our new 
distance estimate of NGC 281, 2.82 kpc, and the so- 
lar motion relative to the LSR based on Hipparcos data 
(Dchncn & Binney 1998), and also assumed the distance 
from the Sun to the Galactic center, R , to be 8.0 kpc 
(Reid 1993), the rotational velocity of the Galaxy at the 
Local Standard of Rest (LSR, at i?o), ©o, to be 236 km s _1 
(Reid & Brunthaler 2004), and a flat rotation curve of 
the Galaxy (i.e. a differential rotation with a nearly con- 
stant rotational velocity O w 0o for the outer Galaxy). 
Subtracting these two effects of the Galactic rotation and 
the solar motion on the apparent motion of NGC 281 
from the observed proper motions, we obtain the mod- 
ified mean absolute proper motions of the maser fea- 
tures in C3 and CI with respect to the Galactic rota- 
tion to be: (nx, A*y)c3-GR =(0.26, —0.62) mas yr _1 and 
(fix, A*r)ci-GR =(-0.98, -2.02) mas yr" 1 , respectively. 
The mean of these two spatial components is expected to 
trace the systemic motion of the NGC 281 region and ob- 
tained to be: (fxx, /Uy)sys-GR =(—0.36, —1.32) mas yr -1 
The direction of this systemic motion with respect to a 
frame rotating with the Galactic rotation is shown as 
a red arrow in figure 3a. At a distance of 2.82 kpc, 
a proper motion of 1 mas yr _1 corresponds to a trans- 
verse velocity of 13.4 km s _1 , and therefore our veloc- 
ity estimate for the systemic motion of NGC 281 is: 
( Vh u 6 )=(-4.8, -17.7)~(-5, -18) km s" 1 in the direc- 
tions toward increasing I and b, respectively (the minus 
signs here thus indicate the motions toward decreasing I 
and b). 

As mentioned in Sato et al. (2007), if we adopt 
the IAU standard values for the Galactic rotation of 
Ro =8.5 kpc and 6 =220 km s" 1 (Kerr & Lynden- 
Bell 1986) instead of R a =8.0 kpc and 9 =236 km s~\ 
and for the solar motion relative to the LSR of 
the velocity V =19.5 km s _1 toward a2ooo =271. °0, 
^2000 =29. °0, then the resulting values vary by up to 
~10% in the b direction: (h'qx, Mgy)= 

(-2.92 -0.25) 

mas yr _1 , (Msx> Msy) = (l-25, —0.61) mas yr _1 and 
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Fig. 3. The NGC 281 supcrbubble on the reproduced diagrams originally by Megeath et al. (2003) (inverse grayscale Hi data from 
Hartmann & Burton 1973; white contour 12 CO data from Dame, Hartmann, & Thaddeus 2001). The positions of NGC 281 (red) 
and IRAS 00420+5530 (green) are indicated, (a) The b vs. I map of the region, velocity-integrated for the Perseus-arm line-of-sight 
velocity range of Vlsr = — 60 to —25 km s — . The measured systemic motion of NGC 281 relative to the rotating Galactic disk 
is plotted as a red arrow, (b) (c) The Vlsr vs - ' diagram, latitude-integrated for the galactic latitude range indicated by a cyan 
rectangle in (a). Two possible rings of clouds with their centers are indicated in cyan and pink (see text). 



((jf x , /U^) sys _GR=(-l-12, -1.28) masyr 1 , yielding a sys- 
temic velocity estimate of (v[ : u£)~(— 15, —17) km s _1 . 

In the I direction, we did not discuss in Sato et al. 
(2007) the velocity component vi parallel to the Galactic 
plane due to its larger dependence of this direction on 
the Galactic rotation model as can be seen in the two 
cases above. The velocity component vi toward decreasing 
galactic longitude obtained above differs from what one 
would expect from a simple expansion of the NGC 281 
superbubble if the origin (e.g., supernova explosions) of 
the superbubble, expected to be around I ~ 121° from 
the ring center in the Vlsr vs. I diagram of figures 3b 
and 3c (see next section), was rotating with the Galaxy. 
The Galactic rotation velocity that would yield vi > (to- 
ward increasing galactic longitude) is 8 > 262 km s -1 for 
Rq = 8.0 kpc; thus the negative velocity component vi 
could be explained by a faster rotational velocity of the 
Galaxy. However, it is likely that this motion vi of the 
region relative to the Galactic rotation is real, reflecting 



the motion of the superbubble as a whole in a large-scale 
phenomenon such as velocity jumps over a spiral shock in 
the Perseus arm of the Galaxy (e.g., Xu et al. 2006; see 
§4.3). 

In the b direction, we find the NGC 281 region at a dis- 
tance of z ~ 311 pc (at d = 2.82 kpc) from the Galactic 
plane and moving away from the Galactic plane with a ve- 
locity component of ~ 18 km s _1 . Assuming the origin of 
the NGC 281 superbubble to lie within the Galactic plane 
(at 6 = 0°), the dynamical timescale of the superbubble es- 
timated from a simple ballistic motion of Vb is: = z/vb = 
bj fiy-GR = (-6.°31)/(-1.36) mas yr^ 1 = 16.7 ~ 17 Myr, 
where /jly—gr lS the derived Y motion of the maser source 
in NGC 281 West relative to the Galactic rotation. This 
timescale, 17 Myr, is consistent with our timescale esti- 
mate of the superbubble expansion in the direction (al- 
most) parallel to the Galactic plane to be discussed in the 
next section. 

Here we also simply estimate the kinetic energy of the 



10 



M. Sato et al. 



NGC 281 region in the direction perpendicular to the 
Galactic plane, and in the next section we will estimate 
the expansion velocity, timcscale, and energy parallel to 
the plane as well. According to the mass estimates by 
Megeath et al. (2002) of the clouds in the region (indicated 
by a cyan rectangle in figure 3a), the velocity-integrated 
mass in atomic and molecular gas are 3.5 x 10 5 M Q and 
10 5 M Q , respectively. Megeath et al. (2002) estimated 
the total kinetic energy of 4.5 x 10 51 ergs for the ring 
of clouds indicated by a cyan (dashed) ellipse in the 
(Z, Vlsr) diagram of figure 3b, expanding at a velocity of 
fring =22 km s _1 . Using the same mass estimate as given 
by Megeath et al. (2002) and our new estimate of the ve- 
locity component i>& ~ 18 km s _1 away from the Galactic 
plane, we estimate the region to also have kinetic energy 
of 3.0 x 10 51 ergs in the direction perpendicular to the 
plane, which requires the energy of multiple supcrnovac. 

Jf.,3. NGC 281 Superbubble: 3D Structure and Expansion 

Our determination of the parallactic distance to 
NGC 281 enables us to discuss the 3D structure and ex- 
pansion of the NGC 281 superbubble. 

In the NGC 281 superbubble, there exists another H2O 
maser souce in the star-forming region IRAS 00420+5530 
("2000 =00 h 45 m , <W> = +55°47' or I =122°.01, b = 
— 7.°07), whose position is indicated by green dots in fig- 
ure 3. Mocllcnbrock, Claussen, and Goss (2007) mea- 
sured the parallactic distance of the H 2 maser source 
in IRAS 00420+5530 with the Very Long Baseline Array 
(VLBA) to be al = 2.17 + 0.05 kpc. Combining the two 
parallax results, one can examine the size, structure and 
the expansion of the superbubble three-dimensionally. 
Figures 4 illustrates the 3D positions of the two sources, 
NGC 281 and IRAS 00420+5530, in (a) edge-on (as a cross 
section at a galactic longitude I ~ 123°) and (b) face-on (as 
seen from the north Galactic pole) views of the Galaxy. 

As mentioned previously, Megeath et al. (2002, 2003) 
found an expanding ring of molecular clouds (called 
'Megeath's ring' hereafter), and estimated from Al (i.e. 
the ring size in galactic longitudes I) to be 270 pc. From 
AVlsr (i- e - the ring 'size' in line-of-sight velocities Vlsr) 
of i> r ingi —22 km s _1 and from A// AVlsr the dynamical 
age is ~ 6 Myr. However, since the heliocentric distances 
to NGC 281 and IRAS 00420+5530 are now determined 
by the parallax measurements to be 2.82 kpc and 2.17 kpc, 
respectively, these clouds (shown as red and green points 
in figure 3) are separated from each by 650 pc (±290 pc) 
along the line of sight. This then implies the major axis 
of Megeath's ring is along the line of sight and a size of 
approximately 700 — 1900 pc in diameter (as derived geo- 
metrically from figure 3b; the radius of Megeath's ring can 
be approximated by the distance between NGC 281 and 
IRAS 00420+5530). However, it remains unclear whether 
such an elongated shape of the molecular ring is plausi- 
ble, compared to the apparent ring size of ~ 300 pc (in 
part of the NGC 281 superbubble with apparent size of 
300 — 500 pc) as seen in the (I, b) plane (figure 3a). 

We thus suggest a smaller ring of clouds indicated by 
a pink (dashed) ellipse (called 'our ring' hereafter) in fig- 



ure 3c, excluding a cloud at Vlsr ^5 ~10 km s _1 from 
Megeath's ring, as another possible interpretation of the 
(Z, Vlsr) diagram. In this interpretation, NGC 281 and 
IRAS 00420+5530 are close to the far side and near side 
of the expanding ring, yielding a smaller ring diameter of 
~ 650 pc along the line of sight, a lower (~70%) expan- 
sion velocity of the ring, T; r i ng 2 — 15 km s , and thus a 
lower (~50%) kinetic energy of 2.1 x 10 51 ergs than for 
Megeath's ring. The sum of this kinetic energy of the ring 
expansion (almost) parallel to the Galactic plane and the 
kinetic energy perpendicular to the plane (estimated to 
be 3.0 x 10 51 ergs in the previous section) yields total ki- 
netic energy of 5.1 x 10 51 ergs, again requiring multiple 
supernovae. We also estimate the timescale of the ring 
expansion parallel to the plane, t r ing2 — 325 pc/i> r ing2 — 
21 Myr, which is in good agreement with the timescale 
tb — 17 Myr independently estimated in the previous sec- 
tion from the velocity component Vb perpendicular to the 
Galactic plane. 

Note that in either interpretation of these two rings, 
the center of the ring (denoted by * in figures 3b and 
3c) is at the same galactic longitude, I ~ 121° (as indi- 
cated by yellow dotted lines in figure 3). Considering the 
ring center as the origin of the superbubble, we can com- 
pare the (observed) line-of-sight velocities of the origin 
(i.e. the ring center in figures 3b and 3c), NGC 281 and 
IRAS 00420+5530 with those expected from the Galactic 
rotation model, expressed by 

Vlsr = #0 sin Zcosfc^-^^ +14 sink, (1) 

where R and © are the galactocentric distance and rota- 
tional velocity of the Galaxy at the source position (I, b) 
in the galactic coordinates, and those for the LSR are Rq 
and Oo Here we also included V z , the source velocity 
component perpendicular to the Galactic plane (toward 
the north Galactic pole) , in order to take into account the 
projection effects on the line-of-sight velocities for sources 
at non-zero galactic latitudes. However, we neglect this 
projection-effect term because sin b ~ —0.1 for NGC 281 
and IRAS 00420+5530 and also because our measurement 
gives V z ~Vb — —18 km s _1 of NGC 281 smaller than the 
galactic rotation (0 and ©0) by an order of magnitude. 
Using the heliocentric distance to the source, d, its projec- 
tion onto the Galactic plane, D' = dcosb, and the cosine 
formula, R is expressed by 

R=^D' 2 + Rl- 2D'R„ cos I (2) 

= yjd 2 cos 2 b + Rl - 2R dcoslcosb. (3) 

Again we adopt 6 ~ ©0 — 236 km s _1 , R ~ 8.0 kpc, 
and then the line-of-sight velocities expected from the 
Galactic rotation model are: V^si = — 36 km s _1 for 
NGC 281, Vy 00420 = -28 km s" 1 for IRAS 00420+5530, 
and Vri ngl = —36 km s _1 for Megeath's ring center (ex- 
pected to be at the same heliocentric distance d— 2.82 kpc 
as NGC 281 from figure 3b), and Vl ng2 = -33 km s" 1 
for our ring center, assumed to be at d = 2.5 kpc. 
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Fig. 4. Schematic diagram of 3D structure and motion of the NGC 281 superbubble: (a) edge-on view (at a galactic longitude of 
I ~ 123°) and (b) face-on view of the Galactic disk. Note that all motion vectors are as seen from the rotating Galactic disk (i.e. 
the Galactic rotation is subtracted). The expected motions are plotted as arrows for the explosion origin (blue dotted; assumed to 
be in the disk), expansion of the superbubble (green dashed) and for the observed motion of NGC 281 (red full). Errors are also 
indicated for our distance estimate of NGC 281 (see text). 



The observed line-of-sight velocities are: Vlsr,N281 = 
-31 km s- 1 for NGC 281 (in 12 CO lines; Lee & Jung 
2003), which is redshifted from the model by AVlsr,N281 — 
+5 km s _1 , and for IRAS 00420+5530, on the other 
hand, Vlsr, 100420 = —51 km s _1 (in C 34 S(3— 2) and other 
lines; Brand ct al. 2001), which is blueshifted from the 
model by AVlsr,ioo420 — —23 km s _1 . This large blueshift 
of IRAS 00420+5530 is likely due to the expansion of 
the molecular ring in the superbubble. The redshift of 
NGC 281 is capable of two different interpretations using 
the Megeath's ring model and ours. 

For Megeath's ring, the line-of-sight velocity of the ring 
center (in figure 3b), VLsR.ringi = —30 km s _1 , is also red- 
shifted from the model by AVLSR,rin g i = +6 km s _1 , sug- 
gesting that not only NGC 281 but also the molecular ring 
as a whole has a redshifted line-of-sight velocity compo- 
nent relative to the Galactic rotation. Then the ring ex- 
pansion velocity is V ring i ~ | AVl S R,I00420 - AVLsR.ringi | - 
29 km s" 1 . 



For our ring, on the other hand, the line-of-sight velocity 
of the ring center (in figure 3c), VLSR,rin g 2 = —36 km s , 
is blueshifted from the model velocity by AVL,SR,ring2 — 
—3 km s _1 , indicating the redshift of NGC 281 is due 
to the ring expansion. Then the ring expansion velocity 
is obtained to be v ring2 ~ |AVLsr,ioo42o - AVL SRiring2 1 = 
20 km s _1 , and the timescale of expansion is roughly: 
iring2 — 325 pc/w r i, lg 2 — 16 Myr, in good agreement with 
our independent timescale estimates above for ring ex- 
pansion along the Galactic disk (i r m g 2 ~ 21 Myr) and for 
vertical expansion out of the disk (tb ~ 17 Myr). 

As mentioned in the previous section, the superbubble 
as a whole seems to have small systemic velocity devia- 
tions (~ 10 km s _1 ) from the Galactic rotation, not only 
along the line of sight (toward blueshift) but also in the di- 
rection toward decreasing galactic longitude. This motion 
relative to the Galactic rotation model is plotted as blue 
dotted arrows in figures 4a and 4b with expected expand- 
ing motion as green dashed arrows and the resulting ob- 
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served velocity of NGC 281 as red arrows. The deviation 
of the superbubble motion from the Galactic rotation, in 
particular the motion toward decreasing galactic longitude 
(against the expected expanding motion toward increas- 
ing longitude), might be reflecting the peculiar motion of 
the Perseus spiral arm, which may be explained by veloc- 
ity jumps over a spiral shock in the Perseus arm. Peculiar 
motion in the Perseus arm is also reported by Xu et al. 
(2006) for the massive star- forming region W3(OH) (the 
position indicated in figure 4b) to be rotating slower than 
the Galactic rotation by ~ 14 km s^ 1 . For the NGC 281 
superbubble, this lends support to the likely peculiar mo- 
tion of the whole superbubble inward the Galactic rotation 
and toward the Sun. 

Also in figure 4b, we have plotted star-forming regions 
in the solar neighborhood and in the Perseus arm of the 
Galaxy with recent parallactic distances precisely mea- 
sured with VERA and the VLBA: Orion KL (Hirota et al. 
2007; Menten et al. 2007; Kim et al. 2008), NGC 1333 
(Hirota et al. 2008), IRAS 16293-2422 in p Oph East (Imai 
et al. 2007) and W3(OH) (Xu et al. 2006; Hachisuka et al. 
2006), demonstrating the high capability of the VLBI 
techniques to reliably trace out the spiral arms of the 
Galaxy. Our parallactic distance to NGC 281 is likely 
to trace the far side of the Perseus arm, and we expect 
to reveal the detailed structure and motion of the Perseus 
arm with future parallax measurements with VERA. 

The parallactic distances to NGC 281 and 
IRAS 00420+5530 revealed the size of the super- 
bubble to be above ~ 650 pc in diameter along the line 
of sight, almost parallel to the Galactic plane. Compared 
to the superbubble size z perpendicular to the Galactic 
plane, z ~ 300 — 400 pc, the superbubble has a structure 
likely elongated along the plane, or at least spherical (i.e. 
by taking the lower limit of 360 pc in size along the line 
of sight) and not elongated vertically in the z direction. 
This is contrary to expectation, i.e. one would expect a 
superbubble to expand and blowout in the z direction 
with a vertical density gradient and with a long timcscalc 
of ~ 20 Myr (see e.g., Tomisaka & Ikcuchi 1986 and 
Mac Low et al. 1989 for superbubble models with vertical 
density stratifaction) . 

It is thus indicated that the expansion of the NGC 281 
superbubble is confined to the Galactic disk. This may be 
explained by Tomisaka (1998)'s 3-dimcnsional MHD sim- 
ulations which show that the superbubble expansion can 
be confined to the disk by the magnetic field of the disk, 
and our estimated 3D structure of the NGC 281 super- 
bubble shows a very analogous shape, size and timcscalc 
to those of model A (size R z ~ 255 pc perpendicular to 
the plane, timescale i ^22Myr) of Tomisaka (1998), which 
assumes a uniform density of the ISM, n = 0.3 cm~ 3 and 
a uniform magnetic field, B(z) = B(0) = 5/iG, aligned to 
one direction within and along the disk. Therefore, the 
shape of the NGC 281 superbubble sugessts that the su- 
perbubble expansion may be confined to the disk due to 
the magnetic field of the disk. 

Within the disk, the Tomisaka (1998) model predicts 
the elongation of the superbubble to be along the direc- 



tion of the magnetic field of the disk, which could be in 
disagreement with the NGC 281 case if the magnetic field 
lies along the Perseus spiral arm (e.g., Han et al. 2006). 
The density decrease of the Perseus spiral arm might be 
also responsible for the structure formation at the position 
and heliocentric distance of NGC 281, thus yielding a pos- 
sible implication on the width of the Perseus spiral arm. 
The further details of the 3D structure of the NGC 281 
superbubble and its relation to the Perseus arm should be 
revealed by future parallax measurements with VERA of 
multiple mascr sources in the superbubble. 

4-4- Distance to NGC 281 and the H-R Diagram 

Our astrometric determination of the distance to 
NGC 281 also provides important information toward bet- 
ter understanding of early-type (O— B) high-mass (M > 
10M Q ) stars and stellar evolution in combination with 
photometric studies of the OB stars in NGC 281 (e.g., 
Henning et al. 1994; Guetter & Turner 1997). 

In the 1990s, the Hipparcos astrometry satellite mea- 
sured the trigonometric parallaxes of more than 100 000 
stars in the solar neighborhood with a precision of about 
1 mas, including direct distance estimates of 20 853 stars 
with a precision better than 10 percent and of 49 399 
stars (in total) with a precision better than 20 per- 
cent (Pcrryman et al. 1995, 1997). Figure 5 is the H- 
R (Hcrtzsprung-Russell) diagram or the color-absolute 
magnitude diagram based on the Hipparcos data of stars 
(shown as black dots) with parallactic distances deter- 
mined with precisions better than 10%. 

The achievements of Hipparcos marked a significant 
milestone in the history of astronomy and brought a new 
era of astrometry, since before Hipparcos, parallax mea- 
surements were limited to the nearest stars and thus di- 
rect distance estimates of only several hundred nearby 
stars were known with precisions better than 20 percent. 
Nevertheless, the parallax measurements by Hipparcos 
were restricted only to within a few hundred parsecs from 
the Sun, which is much smaller than the size of the Galaxy 
(e.g., ~15 kpc in disk radius). 

Due to the low local space density of OB stars, there is 
a lack of early-type high-mass stars with absolute magni- 
tudes reliably derived from astrometrically measured dis- 
tances, as seen on the left side (i.e. the high-mass end) 
of figure 5. Alternatively, most of the distances and thus 
absolute magnitudes of OB stars are estimated by a com- 
parison of photometric colors and magnitudes with an 
empirical standard of zero-age main sequence (ZAMS) 
stars in the color-magnitude diagram, derived from other 
well-studied objects such as Hyades (see Anthony- Twarog 
1982). Such ZAMS comparisons assume compositional 
differences in clusters can be ignored, however, the metal- 
icity dependence of ZAMS stars is important when used as 
distance scale and is the cause of difficulty in a direct com- 
parison between theoretical models and empirical ZAMS 
curves (e.g., Anthony- Twarog 1982; Zinnecker & Yorkc 
2007). Therefore, precise and direct determination of dis- 
tances, and thus absolute magnitudes, are desired and es- 
sential for further studies of stellar evolution sensitivity to 



Distance to NGC 281 



13 




Fig. 5. (a) The H-R (color-magnitude) diagram of stars with parallactic distances determined with precisions better than 10% 
(i.e. ov < O.IO71-). Black dots show the Hipparcos data (Perryman et al. 1997), and red dots show the photometric data of ZAMS 
of early-type stars in NGC 281 by Guetter & Turner (1997), plotted with absolute magnitudes derived with our new parallactic 
distance to NGC 281 of 2.82 kpc. Empirical ZAMS curves are also plotted: Turner (1976) in blue; Blaauw (1963) in pink; and 
Schmidt-Kalcr (1982) in green (dashed). Orange dotted (lower) and dashed (upper) curves show ZAMS curves that would yield 
distances of NGC 281 to be 2.2 kpc and 3.5 kpc (see text), (b) Magnified view of the boxed area (i.e. the high-mass end) in (a). 



metalicity. 

As mentioned earlier, extremely high astrometric accu- 
racy can now be achieved with VLBI techniques out to as 
kiloparsec distances with precisions better than 10 percent 
(e.g., Xu et al. 2006; Hachisuka et al. 2006; Menten et al. 
2007). In particular, recent results of parallax measure- 
ments with VERA are reported in Hirota et al. (2007, 
2008a, 2008b), Honma et al. (2007), Nakagawa et al. 
(2008), Kim et al. (2008), and Choi et al. (2008), clearly 
demonstrating VERA's high capability for the Galaxy- 
scale astrometry. 

Using our new distance determination of NGC 281 by 
parallax measurements with VERA, d = 2.82 ± 0.24 kpc, 
to derive absolute magnitudes, we have plotted the pho- 
tometric data of ZAMS of early-type stars in NGC 281 
by Guetter & Turner (1997) in figure 5 as red dots, which 
give improved coverage of the high-mass end of the H- 
R diagram. In figure 5, the empirical ZAMS curves 
are also plotted for reference: Turner (1976) in blue; 
Blaauw (1963) in pink; and Schmidt-Kaler (1982) in green 
(dashed). Orange dotted (lower) and dashed (upper) 



curves in figure 5 show ZAMS curves that would yield 
distances of NGC 281 to be 2.2 kpc and 3.5 kpc (i.e. the 
distance discrepancy of NGC 281 from previous studies), 
made by additionally shifting the ZAMS by Turner (1976), 
from which Guetter & Turner (1997) obtained a photo- 
metric distance to NGC 281 of 2.94 kpc. 

Figure 5 clearly shows the NGC 281 data agree well with 
the empirical ZAMS curves, and that resolving the large 
distance discrepancy of NGC 281 from previous studies 
was important. The agreement of the parallactic distance 
with the photometric distance is important for confirming 
the accuracy of the correction for extinction or reddening 
in photometric studies such as Guetter & Turner (1997) 
and for enhancing the precision of the H-R diagram. 

Our determination of the parallactic distance to 
NGC 281 therefore contributes to the extension of the H-R 
diagram for the study of high-mass stars and, furthermore, 
opens up new possibilities of using the H-R diagram with 
parallax measurements for Galaxy-scale studies, for ex- 
ample, investigating the formation of the Galaxy via the 
difference in the H-R diagram for components such as the 
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bulge and spiral arms. Further contributions are expected 
with the VLBI techniques and with VERA for direct par- 
allax determination far beyond the solar neighborhood to- 
ward a comprehensive understanding of the Galaxy. 

5. Conclusions 

In this paper, we have presented the results of our multi- 
epoch phase-referencing observations with VERA over 18 
months of an H2O maser source in the Galactic star- 
forming region NGC 281 West, associated with a frag- 
menting supcrbubblc ~ 300 pc above the Perseus spiral 
arm. The primary results are summarized as follows: 

1. We detected a trigonometric parallax of 
0.355±0.030 mas, corresponding to a distance 
of 2.82±0.24 kpc to NGC 281. Our parallactic 
distance agrees well with the photometric distance 
of 2.94±0.15 kpc derived by Guetter & Turner 
(1997), allowing for improved study of the absolute 
magnitutudes of high-mass OB stars, and resolved 
the large distance discrepancy of NGC 281 from 
previous photometric and kinematic studies. 

2. We revised the absolute proper motions of the 
H2O maser features measured in Sato et al. (2007), 
and using our new determination of the parallactic 
distance, we derived more precisely the velocity 
component of the NGC 281 region Vb — 18 km s _1 
perpendicular to the Galactic plane. This yields 
a timescale of % ~ 17 Myr and kinectic energy of 
the region to be 3.0 x 10 51 ergs in the direction 
perpendicular to the plane. 

3. We demonstrated the 3D structure and expansion 
of the NGC 281 superbubbles in comparison of our 
parallactic distance to NGC 281 with the parallactic 
distance derived by Moellenbrock et al. (2007) for 
another H 2 maser source, IRAS 00420+5530 in 
the supcrbubblc. Our new parallactic distance 
revealed the structure of the superbubble ~650 
pc in size parallel to the Galactic disk and with a 
shape slightly elongated along the disk or spherical, 
but not vertically elongated in the z direction. 
Therefore, the superbubble expansion may be 
confined to the disk, possibly due to the magnetic 
field of the disk. 

4. We suggested a new possible interpretation of an 
expanding molecular ring parallel to the Galactic 
plane in the longitude-velocity diagram by Megeath 
et al. (2002, 2003). In either interpretation, the 
ring center, i.e. the likely origin of the superbubble 
lies at I ~ 121°. The velocity deviation of the 
supcrbubblc from the Galactic rotation is estimated 
to be ~ 10 km s _1 inward the Galactic rotation 
and toward the Sun from our proper motion mea- 
surements and from the longitude- velocity diagram, 
which might be due to velocity jumps over a spiral 



shock in the Perseus spiral arm. 

5. We also estimated the velocity component and 
timescale of the ring expansion to be t> r i n g2 — 15 — 
20 km s~ x and t r ing2 ~ 16 — 21 Myr, respectively, 
parallel to the Galactic plane, which are in good 
agreement with those independently estimated for 
the direction perpendicular to the Galactic plane, 
vi, ~ 18 km s _1 and ~ 17 Myr. The kinetic en- 
ergy estimate of the region parallel to the plane is 
2.1 x 10 51 ergs, and the total kinetic energy of both 
perpendicular and parallel to the plane is then esti- 
mated to be 5.1 x 10 51 ergs, requiring multiple su- 
pernovae. 
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